Copolymer adsorption kinetics at a selective liquid-liquid interface: 
Scaling theory and computer experiment 
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We consider the adsorption kinetics of a regular block-copolymer of total length N and block 
size M at a selective liquid-liquid interface in the limit of strong localization. We propose a simple 
analytic theory based on scaling considerations which describes the relaxation of the initial coil into 
a flat-shaped layer. The characteristic times for attaining equilibrium values of the gyration radius 
components perpendicular and parallel to the interface are predicted to scale with chain length 
N and block length M as r± cx M 1+2l/ (here v w 0.6 is the Flory exponent) and as r\\ cx iV 2 , 
although initially the rate of coil flattening is expected to decrease with block size as oc M~ . 
Since typically N 3> M for multiblock copolymers, our results suggest that the flattening dynamics 
proceeds faster perpendicular rather than parallel to the interface. We also demonstrate that these 
scaling predictions agree well with the results of extensive Monte Carlo simulations of the localization 
dynamics. 
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The behavior of hydrophobic - polar (HP) copolymers 
at a selective penetrable interface (the interface which di- 
vides two immiscible liquids, like water and oil, each of 
them being favored by one of the two types of monomers) 
is of great importance in the chemical physics of poly- 
mers. For strongly selective interfaces the hydrophobic 
and polar blocks of a copolymer chain try to stay on 
different sides of the interface leading thus to a major re- 
duction of the interfacial tension between the immiscible 
liquids or melts which has important technological appli- 
cations, e.g. for compatilizers, thickeners or emulsifiers. 
Not surprisingly, during the last two decades the problem 
has gained a lot of attention from experiment 0,12, IE El , 
theory B B-H as well as from computer experiment 
[ToL Ull ll2L Il3l Il4| . While in earlier studies attention 
has been mostly focused on diblock copolymers 0, Q due 
to their relatively simple structure, the scientific inter- 
est shifted later to random HP-copolymers at penetrable 
interfaces Until recently though, the prop- 

erties of regular multiblock copolymers, especially with 
emphasis to their dependence on block length M, have 
remained largely unexplored. 

In a recent study jl5j we showed that the equilibrium 
properties (structure, diffusion coefficient, etc.) of a reg- 
ular HP-copolymer at a selective liquid-liquid interface 
are well described in both regimes of weak and strong 
localization by a scaling theory in terms of the total 
copolymer length N and the block size M. In partic- 
ular, we demonstrated that: (i) the crossover selectivity 
decreases with growing block length as Xc ~ M~( 1+v ^ 2 
and the crossover selectivity to the strong localization 
regime vanishes as Xoo ~ M -1 , and (ii) the size of 
the copolymer varies in the weak localization regime as 
R gJ _ ~ M- 1, ( 1+ "V(i-«') and R g]l ~ M ^-u){\+ v )/(\-u) ^ 

and as R gJ _ ~ M", R g \\ ~ M - ^ - ") in the case of 



strong localization (where v% = 3/4 is the Flory expo- 
nent in two dimensions). We shall use these results in 
the present communication in which we suggest a theory 
of the adsorption kinetics of a regular block copolymer at 
a liquid-liquid interface based on dynamical scaling argu- 
ments. To the best of our knowledge so far there have 
been no attempts to treat this problem analytically or by 
means of computer simulation. 

In our dynamical scaling analysis we consider a coarse- 
grained model of a multiblock copolymer consisting of 
N repeat units which is built up from a sequence of bl- 
and P - blocks each of length M. For simplicity one 
may take the interface with negligible thickness as a flat 
plane which separates the two selective immiscible sol- 
vents. The energy gain of each repeat unit is thus — %, 
provided it stays in its preferred solvent, and the sys- 
tem is considered in the strong localization limit where 
X > Xoo oc M" 1 [H H3. We neglect the effect of bulk 
diffusion on adsorption, and place the center of mass of an 
unperturbed coil at time t = at the interface whereby 
the localization field is switched on. The initial coil will 
start then relaxing with time into a flat ("pancake") equi- 
librium configuration in-plane with the interface and the 
kinetics of relaxation will be determined by the sum of 
the various forces acting on the copolymer. 

In order to estimate the driving force of the flattening 
process, one may recall 0, 0] that the effective attrac- 
tive energy e (per diblock) in the strong localization limit 
is e w X-M. In this case the diblock (i.e. a segment con- 
sisting of one H- and one P-block) plays the role of a blob 
and the overall attractive free energy -F a ttr oc ejV cx xN 
where N ~ N/M is the total number of blobs. Thus the 
effective driving force perpendicular to the interface is 
/attr ~ ~XooN/R± where R± denotes the perpendicular- 
component of the radius of gyration. This force is op- 
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posed by a force of confinement due to the deformation 
of the self-avoiding chain into a layer of thickness R±. 
The corresponding free energy of deformation is simply 
estimated as F r n -nt ~ N(b/R^) 1 / 1 ' where b is the Kuhn 
segment size |l7|. For the respective force then one gets 
f^cNb^/R 1 /^ 1 . 

The equation of motion for R± follows from the con- 
dition that the friction force which the chain experiences 
during the motion in the direction perpendicular to the 
interface is balanced by the sum of /^ tr and f^ nf . In the 
case of Rouse dynamics, each chain segment experiences 
independent Stokes friction so that the resulting equation 
of motion has the form 



CoN- 
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R^M^/ivbN 3 ) so that eqs.(|]J and © can be written in 
the dimensionless form 
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The characteristic times for relaxation perpendicular 
and parallel to the interface in eqs. ©-© should then 
scale as t±_ ~ („b 2 M 1+2ly , and ry ~ ( a b 2 N 2 . 

Eq. JSJ can be solved exactly: 



x 2 (t) 



l-F(2v, 1; l + 2v\ x 1/v (t)) 



x 2 (0) 



I - F{2v,l;l + 2v;x 1/v {$)) 



. (7) 



where Co is the friction coefficient per segment. 

During the flattening process, the chain spreads par- 
allel to the interface due to the excluded volume inter- 
action. Within the Flory mean - field arguments, the 
corresponding free energy F cv ~ vN 2 /(R 2 R±) where v is 
the second virial coefficient and Ru is the gyration radius 
component parallel to the interface. The corresponding 
driving force is /it = vN 2 /(R^R±_). This term is coun- 
terbalanced by the elastic force of chain deformation, i.e. 

by /L = -R\\/(bN). 

Taking into account the balance of forces in the parallel 
direction, the equation of motion for Ru takes then the 
following form 



CoN- 



dR\\ 
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Evidently, the excluded volume interactions provide a 
coupling between the relaxation perpendicular and par- 
allel to the interface. Thus Eqs. and (J2J) describe 
the relaxation kinetics of a multiblock copolymer confor- 
mation at a selective liquid-liquid interface. One may 
readily verify that the equilibrium solutions which follow 
from these equations are 



bM u , 



and 



(vb) llA M~ vlA N v *. 



(3) 



(4) 



This coincides with the equilibrium expressions for R± 
and derived earlier from purely scaling consideration 
(lif. The only difference is in the power of the M - de- 
pendence in eq.Q which looks like M^ v ^ instead of 
M~( v *~ v ) in ref. [15j albeit numerically the values of 
both exponents coincide: v/A ~ (y^ — v) ~ 0.15. 

For the full solution of the equations of motion it is con- 
venient to rescale the variables as x = R±/ (bM v ), y = 



where F(a, /3;7; z) is the hypergeometric function and 
x(0) = R±(0) /bM v is the initial value. In the early stage 
of relaxation (i.e. at t <C tl) one has x 3> 1 and the 
solution, eq.Q, with i?j_(0) ~ bN" reduces to 



Rl(t) = R 2 ± (0)- 



CoM 



(8) 



so that the perpendicular collapse of the chain should last 
proportionally to the block length M. In the opposite 
limit of late stage kinetics, t w t± and x > 1, that is, 
close to equilibrium, the relaxation of R g (t) is essentially 
exponential with t± oc M 2,2 . 



R±(t) 
R 
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Moreover, as typically N S> M, one expects that ry 3> 
t±, i.e. the chain coil collapses first in the perpendicular 
direction to its equilibrium value and after that slowly 
extends in the parallel direction. In this case one can use 
x(t) w x cq = 1 in eq.® and derive the resulting solution 
for the parallel component of R g as 



y(t) - 1 = [j/(0) - 1] exp(-t/T, 



(10) 



where the initial value y(0) < 1. As one may see in Fig. 
^ the time evolution of R± and Rn which follows from 
the solution of the system of differential equations (J5J 
- (0, resembles qualitatively rather well the simulation 
data even though one should bear in mind that the time 
axis for the former is in arbitrary units (scaling always 
holds up to a prefactor) whereas in the simulation time is 
measured in MC steps (MCS) per monomer (i.e., after all 
monomers have been allowed to perform a move at ran- 
dom). The relaxation of Ru looks somewhat slower and, 
as our data for longer chains show, this effect becomes 
much more pronounced with growing number of blocks. 

For the Monte Carlo simulations we use an off-lattice 
bead-spring model that has been employed previously for 
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FIG. 1: Relaxation of R± and Ru with time, following from the 
numerical solution of eqs. ©-© for a copolymer of length N = 
128 and block size M — 16. Both components are normalized by 
their respective equilibrium values. The inset shows the result 
for the same quantities (and the same set of parameters) as 
obtained from the MC simulations. 



simulations of polymers both in the bulk and near con- 
fining surfaces 18]. Recently it was applied by us to 
the study of the static properties of block copolymers at 
a selective interface |l5j, therefore we describe here the 
salient features only. 

Each polymer chain contains N effective monomers 
connected by anharmonic springs described by the 
finitely extendible nonlinear elastic (FENE) potential, 

Ufene = -fR 2 In [1 - ^f^] ■ Here I is the length of 
an effective bond, which can vary in between t m i n < £ < 
tmax, with £ min = 0.4, l max = 1 being the unit of length, 
and has the equilibrium value Iq = 0.7, while R = l ma x — 
£a = £o — imin = 0.3, and the spring constant K is taken 
as K/ksT = 40. The nonbonded interactions between 
the effective monomers are described by the Morse poten- 
tial, U M = £A/{exp[-2a(r-r mm )]-2exp[-a(r-r min )]} 
where r is the distance between the beads, and the pa- 
rameters are chosen as r min = 0.8, em — 1, and a = 24. 
Owing to the large value of the latter constant, Um(t) 
decays to zero very rapidly for r > r min , and is com- 
pletely negligible for distances larger than unity. This 
choice of parameters is useful from a computational point 
of view, since it allows the use of a very efficient link-cell 
algorithm 18]. From a physical point of view, the in- 
teractions Ufene and Um make sense when one inter- 
prets the effective bond as a Kuhn segment, comprising a 
number of chemical monomers along the chain, and thus 
the length unit £ m ax — 1 corresponds physically rather 
to 1 nm than to the length of a covalent C — C bond 
(which would only be about 1.5 A). Since in the present 
study we are concerned with the localization of a copoly- 
mer at good solvent conditions, in Um{t) we retain the 



repulsive branch of the Morse potential only by setting 
Um(t) =0 for r > r mm and shifting U M (r) up by e M - 
The interface potential is taken simply as a step func- 
tion with amplitude x> Ui n t(n, z) — {— <j(n)x/2, z > 
0; er(n)x/2, z < 0} where the interface plane is fixed 
at z = 0, and cr(n) = ±1 denotes a "spin" variable which 
distinguishes between P- and H- monomers. In study- 
ing the flattening dynamics of these chains, we always 
start with a configuration which has been equilibrated in 
a solvent that is good for both types of monomers. At 
time t = 0, the interface is switched on (so that it goes 
through the center of mass of the chain) and the flatten- 




FIG. 2: Variation of the characteristic time for the initial 
relaxation of R± with block length M for chains of length 
32 < N < 512. The theoretical prediction is shown by a dashed 
line. The inset shows the typical behavior of R±(t) in the case 
of N = 256, and M = 2 (solid line), 4 (dotted line), 8 (dashed 
line) and 16 (dot-dash line). 

ing dynamics is then followed until the chain reaches its 
new equilibrium configuration at the interface. Since we 
are interested in the behavior of these chains in the limit 
of strong localization, we have chosen \ = 2\ c where \c is 
the crossover selectivity as obtained in Ref. |l5j . We use 
periodic boundary conditions in the plane of the interface 
while there are rigid walls in the z - direction, where the 
simulation box extends from z = — 32 to z = 32. Typi- 
cally we studied chains with lengths 32 < N < 512 and 
block lengths 1 < M < N/8 whereby all measurements 
have been averaged over 1024 different and independent 
equilibrated starting configurations. 

From Fig. one may readily verify that the initial 
collapse of R±(t <C t±) with time closely follows the pre- 
dicted rate oc M _1 according to eq. © in the limit of 
N,M ^> 1. Moreover, for sufficiently large block size 
M the iV-dependence of this initial rate disappears so 
that asymptotically this initial perpendicular collapse is 
governed by the length M only. In contrast, during the 
late stages of localization at the interface one observes on 
Fig. |21 the expected, cf. eq. 10, characteristic scaling of 
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ti oc M 1+2l/ which progressively improves as the asymp- 
totic limit is approached. Note that the fits of tj_ are 
taken only after the initially unperturbed coil has suffi- 
ciently relaxed, R±(t)/R±(0) < 1/e. As expected, the 
initial dependence of t± on N is seen to vanish as M and 
N become sufficiently large, in agreement with eq. (JSJ. 
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FIG. 4: Variation of tm with chain length N for blocks of size 
M = 2, 4, 8 and 16. The inset shows the typical behavior of 
R\\(t) in the case of M = 4, and chain lengths N = 32 (solid 
line), 64 (dotted line), and 128 (dashed line). The respecive 
equilibrium values are indicated by horizontal lines. 



FIG. 3: Variation of tj_ with block length M for chains of length 
32 < iV < 512. The slope of the dashed line is ~ 2.2, according 
to eq. @. The inset shows the typical behavior of Rj_(t) in 
the case of N = 256, and M = 2 (solid line), 4 (dotted line), 
8 (dashed line) and 16 (dot-dash line). Arrows denote the time 
interval where rj_ has been determined by regression in the case 
of M = 8. 



Eventually, in Fig. 0]we display the measured scaling of 
the characteristic time for spreading along the interface, 
T|| , with block size M for different chain lengths 32 < 
N < 512. Apart from some scatter of data for too small 
J\f = N/M, one recovers nicely the relationship tii oc N 2 
whereby once again, the initial M-dependence gradually 
diminishes as Af 3> 1. 

In summary, it appears that the simple scaling the- 
ory of copolymer adsorption on a selective liquid-liquid 
interface captures well the most salient features of the 
problem as a comparison with extensive MC simulation 
data demonstrates. In this letter we focus our consid- 
eration on the case of Rouse dynamics whereby hydro- 
dynamic effects are neglected. Therefore the present ap- 
proach should be appropriate to the case of an interface 
between immiscible polymer melts. The case of the so 
called Zimm dynamics which accounts for hydrodynamic 
effects will be discussed in an extended paper. 
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